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1.  Establish  mathemat  i  cal  properties  o-f  optimal  in-finite 
dimensional  compensators  based  on  1 i near-quadrat i c- 
gaussi  an  (LG’G)  theory  -for  distributed  systems,  and 
characterise  the  resulting  closed-loop  systems.  Such 
properties  include  the  -form  o-f  control  and  estimator  gains 
and  the  stability  o-f  the  closed-loop  systems.  E:;tend 
previous  results  on  in-finite  dimensional  Riccati  operator 
equations,  which  determine  the  gains  -for  the  in-finite 
dimensional  compensators.  This  work  includes  both 

cont  i  nuoLis-t  i  me  and  discrete-time  control  o-f  -fle:;ible 
structures . 

2.  Develop  approximation  theory  -for  numerical  solution  o-f 
the  in-finite  dimerisional  Riccati  equations  and  construction 
o-f  the  corresponding  compensators.  This  research  concerns 
convergence  analysis  to  determine  necessarv'  and  su-f-ficient 
conditions  for  convergence  of  appro;;imate  solutions  to 
infinite  dimensaional  Riccati  equations  and  the  stu-"1y  of  th 
types  of  convergence  produced  by  different  approximation 
schemes. 

3.  Develop  widely  applicable  formulas  for  computing 
approximations  to  functional  control  and  estimator  gains 


from  finite  dimensional  apprimations  to  the  solutions  to 


in-Finite  dimensional  Riccati  equatior.s. 

4.  Investigate  in-finite  dimensional  autoregressive-moving-- 
average  (ARMA)  models  o-f  distributed  systems  with 
applications  to  large  -fleijible  structures.  This  discrete¬ 
time  representation  o-f  the  distributed  model  o-f  a  -flexible 
structure  will  serve  as  the  basis  -for  adaptive  parameter 

i  dent  1 -f  i  cat  1  on  and  control  o-f  large  space  structures. 

Several  di-f-ficult  existence  and  uniqueness  questions  arise 
in  the  in-finite  dimensional  case  that  do  not  arise  in  -finite 
di men si ons. 

5.  L'evelop  appr  ox  i  mat  i  on  theory  -for  computing  in-finite 
dimensional  ARMA  representations  -for  -flexible  structures 
using  -finite  element  approximations  o-f  the  structures. 


1 1 .  Status  Q-f  Research 

Objectit'e  2:  For  the  continuous-time  case  with  bounded 
input  operator,  the  results  are  parctically  complete. 

Their  application  was  discussed  in  C53  and  will  be 
illustrated  in  L31.  Some  preliminary  discrete-time  results 
are  given  in  LI, 21.  A  compresensi ve  treatment  will  be 
presented  in  L41. 

Objectii'es  2  arid  3:  F'reliminary  results  were  given  in 
C5,61,  and  coiiiplete  re^sults  will  be  contained  in  C3,41. 

Dtjectn-e  4:  Results  are  reported  in  Cl, 21.  A  journal 
paper  extending  these  results  is  planned  -for  SIAM  J.  Contr. 
□pt . 

Objective  b:  While  initial  results  based  on  modal 
appr  o:;  1  mat  i  ori  were  used  in  Cl, 21,  significant  wort  remains 
to  develop  theory  and  methods  for  arbitrary  finite  element 
a  p  p  f  D ; .  i  m  a  1 1  o  n  s; . 

Fcest^arch  in  discrete-time  control  and  parameter 
estimation  of  flexible  space  structures  is  continuing  under 
AFOGFc  Grant  840309,  "Optimal  Control  and  Identification  of 
Large  Space  Structures:  Distributed  System  Theory  and 
Numeircal  Approximation."  Also,  prpblems  with  unbounded 
input  operators  are  being  investigated,  especially  in  the 
discrete-time  case. 
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Summary 

Approxlmatior.  of  an  ideal  infinite  dimen¬ 
sional  compensator  for  an  infinite  dimensional 
system  is  discussed.  Tvo  approaches  are  con- 
slderedi  one  is  to  design  a  compensator  for  a 
finite  dimensional  model  of  the  system,  and  the 
other  approach,  which  is  stressed  here.  Is  to 
truncate  an  infinite  dimensional  ARMA  representa¬ 
tion  of  the  ideal  compensator.  Performance  of 
the  closed-loop  system  is  discussed,  and  an 
example  from  control  of  flexible  structures  is 
presented . 


Introduction 

This  paper  treats  compensator  design  for 
a  discrete-time  infinite-dimensional  linear 
control  system.  Usually,  the  compensator  which 
solves  the  optimal  LQC  problem  for  such  a 
system  is  infinite  dimensional;  in  other  words, 
it  has  an  Irrational  transfer  function.  We  will 
refer  to  an  infinite  dimensional  compensator  as 
an  ideal  compensator  because,  in  applications, 
only  finite  dimensional  compensators  can  be 
built,  preferably  of  low  order.  The  design 
philosophy  of  this  paper  is  to  construct  finite 
dimensional  compensators  which  approximate  the 
ideal  infinite  dimensional  compensators. 

For  the  discrete-time  problem  discussed 
bere,  we  present  two  approaches  to  approximating 
an  infinite  dimensional  compensator.  The  first 
approach  is  to  approximate  the  infinite  di¬ 
mensional  control  system  by  a  finite  dimensional 
model  and  then  design  a  corresponding  finite  di¬ 
mensional  compensator.  We  call  the  resulting 
implementable  compensator  an  indirect  approxima¬ 
tion  of  the  ideal  compensator.  This  design 
procedure  is  common  practice  among  engineers 
who  use  finite  dimensional  control  methods  for 
flexible  struetures-and  often  works  quite  well 
-although  the  finite  dimensional  compensator 
Is  rarely  thought  of  as  approximating  some 
ideal  compensator.  Our  only  suggestion  is  to 
try  to  ensure  that  the  finite  dimensional 
compensator  will  perform  well  by  selecting  the 
approximation  scheme  so  that,  as  the  dimension 
of  the  model  increases,  the  finite  dimensional 
compensator  converges,  in  an  appropriate  sense, 
to  the  ideal  Infinite  dimensional  compensator. 

The  second  approach  to  designing  s  finite 


dimensional  compensator  produces  what  we  call 
a  direct  approximation  of  the  ideal  compen¬ 
sator.  In  this  approach,  an  infinite  di¬ 
mensional  ARMA  representation  of  the  ideal 
compensator  is  constructed  and  then  truncated 
to  obtain  the  ARMA  representation  of  the  finite 
dimensional  compensator.  For  computational 
purposes,  the  ARMA  representation  of  the 
infinite  dimensional  compensator  is  obtained 
as  the  limit  of  the  ARMA  representations  of 
the  finite  dimensional  compensators  obtained 
by  the  indirect  approach.  But  these  compen¬ 
sators  are  discarded  once  the  ARMA  coefficients 
for  the  infinite  dimensional  compensator  are 
obtained,  and  the  ideal  compensator  itself  is 
truncated  to  obtain  the  direct  approximation. 

In  this  paper,  we  are  most  interested 
in  the  direct  approximation  of  the  ideal 
compensator.  We  will  discuss  the  performance 
of  the  closed-loop  system  with  the  resulting 
finite  dimensional  compensator,  showing  that 
stability  and  near-optimality  can  be  preserved. 
Finally,  we  give  an  example  in  which  a  flexible 
beam  is  controlled. 


The  Infinite  Dimensional  Compensator 

Let  E  be  an  infinite  dimensional  Hilbert 
space,  T  a  bounded  linear  operator  on  E,  B  a 
linear  operator  from  R*  to  E,  C  a  bounded 
linear  operator^from  E  to  R**,  and  D  a  linear 
operator  from  r"  '.o  R^.  We  consider  the 
discrete-time  control  system 

x(t*l)  ■  T  s(t)  ♦  B  u(t),  f0,l,2 .  (1) 

y(t)  ■  C  s(t)  ♦  D  u(t)  ,  (2) 

where  the  state  vector  s(t)eE,  the  control 
vector  u(t)t  r"  and  the  measurement  vector 
y(t)E  rP. 

Our  ideal  compensator  for  this  system  is 

on  the  infinite  diMennionnl  di screte^tiae 
•tnte  estiaator 

£{tel)  ■  T  8(t)  ♦  Bu(t)e  C(y(t)  -C8(t)  -  D  u(t)),  (3) 

where  the  estimator  gain  G  is  a  linear  operator 
from  R**  to  E.  The  feedback  control  is 

u(t)  -  -Ks(t)  (I.) 


This  research  was  supported  by  HSF  Grant  where  the  control  gain  K  is  a  bounded  linear 
EHGTS-ObTJl.  operator  from  E  to  R  .  Our  first  representa¬ 

tion  of  the  ideal  eomnensator  is  aiven  by  13) 


>na  ( li ) .  We  vill  call  this  compensator  optimal 
If  C  and  K  arc  chosen  accordlns  to  an  infinite- 
dimensional  dlsereto-time  LQO  problem  for  (l) 
and  (2)  (with  added  noise).  This  Is  not  nee- 
cesary  for  the  discussion  In  this  paper,  al¬ 
though  we  have  chosen  C  and  K  this  way  In  the 
Example.  We  call  the  compensator  in  (3)  -  (^) 
Ideal  because  It  Involves  Infinite-dimensional 
dynamics. 

'We  will  assume 


The  Closed-loon  System  with  the  ARMA  Compensator 

Henceforth,  ve  will  assume  that  the 
sa^uenees  u{-j)  and  y(-j),  J  ■  0,1,2,...,  are 
In  t-(R  )  and  respectively.  In  other 

words,  the  control  and  measurement  histories 
arc  cquarc-summable .  Thus,  If  z(t)  is  given 
by  (10),  (6)  holds.  For  t  b  0,  ve  define 

U(t)  «  {u(t),  u(t-l),  tt{t-2),...)t  t2(R“)  (ll.) 


Hypothesis  1.  The  spectral  radius  of  T-CC  Is 
less  than  1.  □ 

This  means  that  there  exist  constants  M  and  r, 
r  <  1,  such  that 

II  (T-CC)"1|  s  Mr",  n  •  1,2,...  .  (5) 

In  applications  to  flexible  structures. 
Hypothesis  1  implies  that  all  but  a  finite 
number  of  modes  are  uniformly  exponentially 
stable  due  to  inherent  system  damping,  and  that 
any  unstable  modes  are  observable.  Also,  from 
(  ^  )  it  follow  r,  that  t(t)  approaches  2(t)  expo¬ 
nentially;  i.c., 

II  t(t)-z(t)l|  S  Mr^ll  z(0)-2{0)l|  ,  t-1,2,...  . 

(6) 

To  derive  the  infinite  dimensional  ARMA 
representation  of  (3)  -  -  ( )  ,  ve  define 

S  -  T  -  CC  and  B  -  B  -  GD,  (7) 

and  write  (3)  as 

z{tel)  ■  S  z(t)  ♦  B  u(t)  ♦  C  y(t).  (8) 

Repeated  application  of  (8)  yields 

z(t)  -  t  S^"^(B  u(t-j)  ♦  0  y(t-j))  ♦  s"t{t-n), 

(9) 

n  w  1,2,...  . 

If  we  set  z{t)  ■  0,  u(t)  «  0,  and  y(t)  »  0  for 
t  s  0,  then  from  (9)  ve  have 

z(t)  »  t  u{t-j)  ♦  c  y(t-j)),  0  <  t< 


T(t)  ■  (y(t),  y(t-l),  y(t-2),...}c  tgCR**)  (15) 

It  is  now  an  easy  exercise  to  show  that,  when 
the  compensator  in  (11)  is  used,  there  exists 
a  constant  such  that 

II  -K  z(t)  -  u(t)||  S  M^r^dl  i(0)||  ♦  II  U(0)||  ♦  ||  Y(0)||  )  , 

t  •  1,2,...  ,  (16) 


II  (T-BK)‘z(0)  -  z(0)||  S  Mjr‘(||  t(0)||  +  ||  U(0)||*  ||  Y(0)||), 

t  -  1,2 .  (17) 

For  direct  approximation  of  the  ideal 
compensator,  it  will  be  useful  to  define  the 
space  Z  «  E  X  t2(R")x  and  to  define  the 

unique  bounded  linear  operator  T  on  Z  such 
that,  for  t  ■  0,1,2 . 

(z(t+l),U  (t+1),  Y(t+1))  -  T(z(t),  U(t),  Y(t)).  (18) 

It  would  be  easy  to  write  out  T  explicitly, 
but  this  should  not  be  necessary.  Let  us  note 
only  that  if,  in  addition  to  Hypothesis  1, 

T-BK  is  uniformly  exponentially  stable,  then 
so  is  T;  i.e.,  the  spectral  radius  of  T  is 
less  than  1.  The  purpose  of  (16)  is  to  allow 
us  to  view  the  closed-loop  system  as  an 
autonomous  discrete-time  process  evolving  on 
the  state  space  Z.  When  the  ARMA  compensator 
is  used,  the  closed-loop  system  at  time  t  is 
the  vector  (z(t),  U(t),  T(t))E  Z;  i.e.,  the 
vector  consisting  of  the  current  state  of  the 
control  system  in  (l)  and  the  current  histories 
of  the  control  and  the  measurement. 


Fur'.-hermore ,  if  the  sequences  rJ||z(-j)||  , 

r'^||“(-j)  11  .  rdly(-j)II,  (j  ■  0,1,2,^..) 

are  summable,  then  ve  still  can  define  z(t)  by 
(10), and  (6)  will  hold. 

With  (10),  the  control  lav  in  (b)  becomes 

S  OD 

u(t)-  -Z  a  u(t-j)  -r  8,  y(t-j),  t-  0,1,2 .  (li) 

J-1  J  J-1  ^ 

where  each  Oj  is  an  mxm  matrix  and  each  8.  is 
and  mxp  matrix.  The  control  lav  in  (ll)  *is 


Indirect  Approximation  of  the  Ideal  Compensator 


The  idea  here  is  to  approximate  the 
infinite  dimensional  control  system  in  (1)  and 
(2)  by  a  sequence  of  control  systems,  each  of 
finite  dimension  N.  Thus  ve  vill  have  se¬ 
quences  of  approximating  operators  T„ ,  B„,  C„ , 

Iv  N  Iv 

and  D||.  For  each  H,  an  H-dimensional  compen¬ 
sator  of  the  form  (3)  -  (^)  is  designed  with 
estimator  and  control  gain  s  0,  and  K„.  re- 


the  ARMA  representation  of  the  ideal  compensator  spectively.  This  H  compensator  it  then 


in  (3)  - 
by 


L  The  BAtricet  a.  and 
*  J 


are  given 


a.  ■  K  S'*"^B 


knd  8 K  S^'^G, 
J 


and  in  view  of  ($)  and  (7),  ve  have 
II  »  II  <  II  rll .  II  nil  .  u  II  A  II  <  III 


ijll  s  11  Kll*  11  Bll  •  M  r**"^  and  ||  8^11  S  INI*  H  g1|  'Mr^ 

(13) 


Zu(t*l)  ■  T||Zu(t)  ♦  B||U(t)  ♦  C||(y(t)  -  C||Z||(t)-DjjU(t)), 

(19) 

u(t)  -  -KjSj(t)  .  (20) 

If  the  approximation  is  done  correctly 
(see  [1,2]  for  related  continuous-time  problems 


and  [3l  for  the  dlacrete^tine  enae),  then  the 
conpenaator  in  (19)  -  (20)  vlll  anprozlaate 
the  ideal  eeapenaator  In  (3)  -  (^)  an  I  beeoaea 
large.  Although  atronger  conTergenee  often  can 
be  obtained,  here  we  will  need  only 

Hypotheala  2.  Let 


(a(t+l).  U(t^l),  r(tn))_  _ 
-"l*"2 


.  ,  <27) 

■  T  ■  (*(t).  U(t),  7{t)) 

’l*"2 

low  let  ua  write  (ll)  and  (26)  aa 


- 


and  B. 


- 


Sj,  ♦  S  atrongly, 

B||  ♦  B  strongly, 

Ojj  ♦  C  strongly, 

and 

K„  K  strongly.  □ 


(21)  u(t^l)  -  KiU(t)  ♦  KjKt)  (2B) 

and 

u(t+l)  -  K  U(t)  ♦  r.,  Y(t),  (29) 

(22)  1  *2 

(23)  reapeetively ,  where  C  ,f.‘)  and 

(2J.)  Kg*  ''^2B2®  <<^'<2^’’^*  From  (13),  it 

follows  that 

(25)  11  Kj  -  K^g^II  -  0  as  Bj  *  -  (30) 


The  statement  that  the  finite  dimensional 
compensator  in  (19)  -  (20)  converges  to  the 
infinite  dimensional  compensator  in  (3)  -  (^) 
means  that  S_,  B  ,  ,  and  K  converge  in  some 

n  n  Q  O 

sense  to  S,  fi,  C,  and  K,  reapectivily .  Also, 
under  additional  but  often  realistic  hypotheses, 
it  can  be  shown  that  the  closed-loop  system 
consisting  of  (l),  (2),  (19),  and  (20)  converges 
to  the  ideal  closed-loop  system  consisting  of 
(l),  (2),  (3),  and  (U)  in  such  a  way  as  to 
preserve  stability  and  near -optimality  for  > 
sufficiently  large.  However,  in  this  paper 
we  will  pursue  such  convergence  only  for  the 
subsequent  direct  compensator  approximation. 


and 

II  Kg  -  Kgg^ll  -  0  as  gg  *  “•  <31) 

Finally,  it  is  straightforward  to  show 

■>  *  I*  V*^2Ngll  *  <32) 

so  that 

11  T  -  *•  "i-’a  *  •  •  <33) 


Direct  Approximation  of  the  Ideal  Compensator 

Truncation.  The  idea  now  is  to  compute,  as 
discussed  later,  the  matrices  and  6^  for 

the  ARMA  compensator  in  (ll),  and  then  to  ap¬ 
proximate  this  compensator  directly.  The 
obvious  way  to  approximate  (11)  is  simply  to 
truncate  all  but  a  finite  number  of  terns  in 
each  of  the  series  in  (ll).  That  is,  replace 
(ll)  by  the  control  law 


u(t)  -  -I  a  u(t-j)  -  I  6  y(t-j),  t  ■  0,1,2,...  .  (26) 
3-1  ^  J-1  ^ 

For  any  positive  integers  and  Rg,  (26)  is 

the  finite  dimensional  ARMA  approximation  of 
the  ideal  compensator. 

Converxence  and  Performance  of  the  Closed-loop 
System. 

When  the  ARMA  compensator  in  (?6)  is  used 
to  control  the  system  in  (l),  we  again  will 
call  the  triple  (s(t),  U(t),  T(t9_  _  the 

■l*"2 

closed-loop  system  at  time  t.  lots  that  U(t) 
and  T(t)  are  still  the  infinite  histories  in 
(!>))  and  (15),  so  that  (s(t),  U(t),  T(t))is  a 
vector  in  Z  whether  (11)  or  (26)  is  used. 

Also,  we  can  define  a  bounded  linear 
operator  T  on  Z  such  that  the  closed-loop 

"l  ’  2 

system  resulting  from  (26)  can  be  written 


From  (33),  we  see  that  the  response  of 
the  closed-loop  system  obtained  by  using  the 
finite  dimensional  compensator  in  (26)  can  be 
made  arbitrarily  close  to  the  response  of  the 
ideal  closed-loop  system  obtained  by  using 
the  ideal  compensator  in  (ll).  if  the  ideal 
closed-loop  system  is  uniformly  exponentially 
stable,  then  so  is  the  closed-loop  system  with 
the  finite  dimensional  ARMA  compensator  for 
H^  and  Hg  sufficiently  large.  Also,  if  the 

ideal  closed-loop  system  is  optimal  in  some 
sense,  then  the  closed-loop  system  obtained 
with  (26)  can  be  made  arbitrarily  closed  to 
optimal.  Of  course,  we  hope  that  in  applica¬ 
tions,  stability  and  near  optimality  can  be 
obtained  with  reasonably  small  and  Ng. 


Computation  of  a  and  6 

-  ^ 


ll 


For  the  finite  dimensional  state  estimator 
and  control  law  in  (19)  and  (20),  we  can 
retrace  the  development  in  (7)  -  (11)  to 
obtain  an  approximation  of  (11)  with  a  and 

replaced  respectively  by 

<3'*) 

A 

Where  and  B^  are  given  by  (21).  In  other 
words,  after  computing  the  gains  C||  and 

for  (19)  and  (20),  we  replace  the  operators 
K,  0,  8,  and  B  la  (12)  by  Ith  order  approxima¬ 
tions.  Hypothesis  2  then  guarantees 


o,_  ♦  a,  *nd  6,_  ♦  S,  »•  ■  ♦  .  (35) 

J  H  J  J  W  J 


Ad  Infinite  Dlaenslonal  MA  Conpenaator 


In  general,  stability  of  both  T-GC  and 
T-BK  does  not  guarantee  stability  of  the 
operator 

S  -  T  -  CC  -  BK  -  ODK  (36) 

However,  S  often  Is  uniformly  exponentially 
stable.  When  this  Is  the  case.  It  is  conven¬ 
ient  to  write  (3)  and  ( U )  as 

r(t+l)  ■  §  t(t)  ♦  Qy(t),  (37) 

u(t)  -  -K  J(t)  .  (36) 

Then,  retracing  the  steps  to  (ll)  leads  to 


“(t)  -  y(t-3),  t  -  0,1,...,  (39) 

where  each  y.  is  an  mvp  matrix  given  by 
J 

Yj  -  K  S-^'^G  .  (ItO) 

The  control  lav  in  (39)  is  the  MA  representa¬ 
tion  of  the  ideal  compensator. 

or  course,  we  approximate  (39)  by 


u(t)  •  -Z  Y,  y(t-J  )  .  (Ill) 

J-1  J 

The  analysis  of  the  convergence  and  performance 
of  the  closed-loop  system  produced  by  (26)  is 
easily  modified  to  obtain  the  same  results 
when  the  finite  dimensional  MA  compensator  in 
(ill)  is  used.  Also,  the  matrices  Y. 
computed  as  the  limits  of  ^ 


*1-1 


Example 


We  will  control  the  transverse  vibrations 
of  a  simply  supported  Euler-Bernoulli  beam 
represented  by  the  partial  differential  equa¬ 
tion 

w(t,n)  ♦  2CA^^^  w(t,n)  ♦  A„w(t,n)  -  o  (U3) 

o  o 

where  n  c(0,l)  is  measured  along  the  beam. 

The  basic  space  for  v(t,*)  is  L2(0,l),  and 
the  operator  A  is  defined  by 


D(A^)  -  (♦  C  H  (0,1):  ♦(O)  -  ♦(!)  ■  ♦"(O)  -  *"(1)  .  O) . 


A  <(>  - 
o 

We  take  the  damping  ratio  C  >  .01. 


Our  control  ii 


of  the  beam.  Hence  .^he  boundary  condition 


8  w(t.n) 


This  control  is  a  piecewise  constant  function 
of  time  t,  and  the  length  of  the  time  step  is 
.03. 

To  put  the  system  In  first  order  form, 
we  write  s(t)  -  (v,(r)  and  take  the  space  E 

to  be  H^(0,1)  X  L^(0,1),  where  H^(0,1)  - 

(♦  E  4(0)  ■  4(1)  “  0).  Then  (Ii3)  becomes 


where  A  (see  [1])  generates  a  C  -semigroup 
T(*)  on  E.  The  operator  T  in  °(l)  is  T(.03) 
and  the  operator  B  can  be  obtained  from  the 
solution  to  (k3)  for  constant  u(t)  in  (U6). 
Note  that  B  is  bounded  for  the  discrete-time 
problem,  although  the  input  operator  is  un¬ 
bounded  for  the  continuous-time  problem. 


Our  measurement  is  the  slope  at  the  left 
end  of  the  beam: 


8w( t .n ) I 


This  operator  C  is  a  bounded  linear  functional 
2 

on  E  >  H  xLg. 

For  the  estimator  gain  G  we  solved  for 
the  optimal  filter  for  (l)  and  (2)  when  the 
process  noise  for  (1)  has  covsriance  operator 
I  and  the  measurement  noise  for  (2)  has 
covariance  1.  Because  the  resulting  estimator 
had  unsatisfactory  eigenvalues  (the  eigenvalues 
of  S  in  (T)),  ve  added  an  o-shlft  to  obtain  a 
G  which  made  the  eigenvalues  of  S  all  have 
magnitude  less  than  .9. 

For  the  control  gain  K,  we  solved  the 
optimal  regulator  problem  for  (l)  and  the 
performance  index 

J  “klof  II  He  ♦  “^(-o^k))  .  (1.9) 

To  compute  G  and  K  ve  approximated  the 
solutions  to  infinite  dimensional  discrete¬ 
time  Riccati  equations,  as  discussed  in  (3). 
Basically,  ve  projected  the  infinite  dimen¬ 
sional  LQG  problem  on  subspaces  spanned  by  N 

natural  modes  of  the  beam,  H  ■  1,2 .  Tbe 

ARMA  gains  ,  and  Yj  were  computed  as  tbe 

limits  of  (3I1)  and  (U2)  for  increasing  N. 

We  found  that  §  (see  (36))  is  exponentially 
stable  for  this  example. 

Tbe  following  tables  give  the  values  of 
the  ARMA  and  MA  compensator  gains. 


oment  on  the  left  end 


1 

-.019 

-.020 

-.020 

1 

-1.387 

-1.38lt 

-1.383 

2 

.  lilt 

.117 

.119 

2 

-  .6Ii3 

-  .638 

-  .638 

3 

.  0lt5 

.0li9 

.051 

3 

1.  It  It  It 

l.blt2 

l.ltItU 

U 

.027 

.  022 

.021 

It 

.615 

.620 

.623 

5 

.  198 

.193 

.191 

5 

-2.11t7 

-2.151 

-2.15lt 

. 

. 

• 

• 

. 

• 

• 

. 

. 

. 

. 

• 

• 

• 

• 

• 

15 

.  029 

.  028 

.026 

15 

.080 

-  .082 

-  .079 

16 

.  03I1 

.  03lt 

.  03lt 

16 

-  .2llt 

-  .221 

-  .223 

17 

-.02lt 

-.023 

-.023 

17 

.  050 

.0lt9 

.0lt8 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

. 

. 

• 

• 

• 

• 

. 

. 

. 

. 

• 

• 

• 

• 

35 

-.001 

-.001 

-.001 

35 

.  007 

.  008 

.  008 

36 

.  003 

.  003 

.003 

36 

-  .010 

-  .008 

-  .008 

37 

.001 

.  001 

.001 

« 

37 

.  001 

.  0001 

-  .0001 

• 

• 

• 

• 

; 

• 
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Abstract 


This  paper  presents  an  Infinite  dlaen- 
slonal  ARHA  aiodei  for  Infinite  dlnenslonal 
tlBe>1nvarlant  linear  systeais  In  which  the 
state  transition  operator  Is  trace  class. 

The  ARHA  Model  appears  to  be  the  natural 
extension  of  the  Hnite  dlaenslonal  ARHA 
aiodei.  and  It  Is  coiiputed  as  the  Unit  of  a 
sequence  of  finite  dlaenslonal  ARHA  aiodels. 

In  a  flexible  structure  exaaiple.  nusierlcal 
results  are  given  for  tKe  ARHA  aodels  of  both 
the  epen>loop  plant  and  an  optlaial  coapensator. 

I .  Introduction 

In  [4],  we  Introduced  an  Infinite  diaien* 
slonal  ARHA  aodel  for  a  class  of  Infinite 
dlnenslonal  d1screte>t1ns  control  systens. 

The  method  there,  based  on  an  infinite  dimen* 
slonal  observer  for  the  state  apace  represents* 
tion  of  the  system,  can  produce  en  Infinite 
number  of  different  ARHA  models  of  the  sane 
system,  but  none  of  those  ARHA  models  Is  a 
natural  extension  of  the  minimal  order  ARHA 
model  for  a  finite  dimensional  system  repre* 
sented  In  state  space  form.  For  a  large  class 
of  Infinite  dimensional  lystems,  we  present 
In  this  paper  the  natural  extension  of  the 
minimal  order  finite  dimensional  ARHA  model. 
This  paper  deals  exclusively  with  single* 
1nput/t1ngle*output  systems,  but  the  main 
Ideas  and  results  can  be  extended  to  the 
mu1t1*1nput/mu1t1*output  case. 

That  the  Infinite  dimensional  ARHA  model 
of  this  paper  Is  the  natural  extension  of  the 
finite  dimensional  case  should  be  clear  from 
the  way  the  present  ARHA  model  Is  defined  In 
Section  2  and  from  certain  of  Its  properties. 
One  of  the  most  Important  properties  Is  that 
the  ARHA  coefficients  corresponding  to  the 
Infinite  dimensional  state  space  representa* 
tion  of  the  system  are,  under  appropriate 
convergence  criteria,  the  limits  of  the  ARHA 
coefficients  corresponding  to  Increasingly 
large  finite  dimensional  approximation  of  the 
Infinite  dimensional  system.  See  Section  4. 

After  defining  the  ARHA  model  and  giving 
some  of  Its  Important  characteristics  In 
Sections  2  and  3,  we  sketch  the  relevant 
approximation  theory  In  Section  4.  geceute  ef 
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of  limited  space,  we  will  not  prove  all  of  the 
results  stated  here,  but  we  will  Indicate  the 
main  Ideas  of  the  missing  proofs.  In  Section 
S,  we  apply  the  results  to  obtain  the  ARHA 
awdel  ef  a  flexible  structure  and  the  ARHA 
representation  of  an  optimal  (Infinite  dimen* 
slonal)  LQ6  compensator  for  the  structure. 

At  the  example  Illustrates,  we  are  Inter* 
ested  In  both  the  ARHA  model  for  the  open*leop 
version  of  a  distributed  parameter  control 
system  and  the  ARHA  model  for  an  Infinite 
dimensional  compensator  for  auch  a  system. 

The  ARHA  model  of  the  open* loop  system  should 
be  useful  for  Infinite  dimensional  adaptive 
control  and  parameter  Identification  theory, 
and  at  a  limit  that  reveals  tha  limiting 
properties  of  the  ARHA  representations  of 
finite  dimensional  approximations  to  the  dis¬ 
tributed  system.  The  ARHA  representation  of 
an  Infinite  dimensional  compensator  for  a  dis¬ 
tributed  system  gives  a  concrete  representa¬ 
tion  of  an  otherwise  abstract  entity,  and  can 
serve  as  a  basis  for  deriving  Implementable 
(finite  dimensional)  approximations  to  such  a 
compensator. 

He  consider  control  systems  whose  state 
space  representations  have  the  form 

i(trl)  ■  Ti(t)  ♦  $  u(t),  t  •  0,1,2...  (1.1) 

y(t)  •  C  i(t),  (1.2) 

where  the  state  vector  i(t)  Is  In  adlllbert 
space  E,  T  Is  a  bounded  linear  operator  en  E, 

$  Is  a  linear  operator  from  the  real  line  Into 
E  (1.0.,  $  t  E)  and  C  Is  a  bounded  linear 
functional  on  E.  The  scalar  H(t)  Is  the 
control,  and  the  scalar  y(t)  Is  the  measurement. 
To  define  the  ARHA  model  of  this  paper,  we 
must  assume  that  the  operator  T  Is  trace  class 
(see  [2,S]  and  Section  2.  Also,  we  assume  that 
the  system  (1.1)*(1.2)  Is  observable  In  the 
sense  that  there  Is  no  nentere  s(0)  such  that 
y(t)  ■  0  for  each  t  a  0. 

The  condition  that  T  be  trace  class  holds 
for  the  discrete-time  state  space  representa¬ 
tion  of  the  solutions  to  many  general  lied  wave 
equations,  heat  oquatlens  and  time-delay 
equations.  As  discussed  In  Section  S,  a  common 
modal  damping  model  ensures  that  T  Is  trace 
class  for  the  linear  distributed  model  of  a 
flexible  structure. 


Undtr  the  above  hypotheeei.  the  ARNA  aodel 
of  the  lyttea  (1.1)>(t-2)  hat  the  fora 


(1.3) 


y(t)  •  I  (e,y(t-l)  ♦  bfU(t-1)j. 
l-l 

■here  the  ARHA  coefficlentt  a^  aitd  h^  are 

tellers.  The  purpose  of  this  paper  It  to 
define  these  coefficients  so  that  the  aeasure* 
■ent  In  (1.2)  does  satisfy  (1.3).  to  establish 
the  decay  rate  for  e^  end  b^  alth  Increasing  1, 

to  Indicate  hov  the  coefficients  can  be  tea- 
puted,  and  to  illustrate  the  results  vlth  a 
nuaerical  exaaple. 


y(t)  In  (1.2)  Indeed  satisfies  (1.3)  requires 
a  nontrivial  proof  along  with  precise  definl* 
ttons  of  the  Infinite  output  and  Input 
histories  essuaed  In  (1.3).  leceuse  of  Halted 
•pace.  «e  only  «111  outline  this  proof  In 
Section  4. 


111.  Decay  Rates  for  the  ARWA  Coefficients 


Theorea  1 .  For  each  r  >  0,  11a  r^*.  •  0. 

1«o  ' 

In  other  aords.  the  a^'s  decay  fatter  than  any 
oaponentlal. 


II .  Definitions  of  the  ARHA  Coefficients 


2.1.  Trace  Class  Operators 


The  proof  of  Theorea  1  requires  the 
folloving  leaaa.  •hich  It  the  generalization 
of  the  finite  diaentlonal  situation  and  can 
be  proved  using  (2.2). 


He  vlll  need  toae  standard  results  on 
trace  dess  operators  on  Hilbert  spaces. 

(See  [2.  p.1088-n0S]  and  [5.  p.S21-S24].) 

First,  a  linear  operator  Is  trace  class  If  It 
It  coapact  and  Its  singular  values  are  tuaaable. 
In  this  case,  the  trace  nora  of  T,  denoted  by 


Leaaa  1 .  For  1  ■  I.2.***.  the  coefficient 
a^  is  (-1)^*'  tiaes  the  absolutely  convergent 

countable  sua  of  all  products  of  1  distinct 
eigenvalues  of  T. 


.  Is  equal  to  the  sua  of  the  singular 


values  repeated  according  to  aultipllcity. 
Froa  here  on.  T  vlll  be  trace  class. 


Proof  of  Theorea  1.  For  1  •  1.2.»*»,  let  the 

ets 


products  of  1  distinct  eigenvalues  of  T  be 


Let 


ordered  and  denoted  by  k  ■  I.2.***.  so 

that 

»(1) 


of  T.  Then 


tr  T  -  I  L, 


(2.1) 


a,  -  I  n'’'  .  1  •  I.2.- 

1  Ip)  » 


(3.1) 


1-1 


Peat,  define 


and 


•  I  •  »  •  I.*.*”. 


Ilk, I  *  l|T||,  < 
1«1 


(2.2) 


(3.2) 


The  set  of  trace  class  operators  on  a  Hilbert 
space  E  Is  a  linear  space,  and  the  coaposltlen 
of  a  trace  class  operator  «1th  a  bounded  linear 
operator  Is  trace  dess. 


that  each  Hs  finite  can  be  seen  by  defining 

the  Infinite  diaentlonal  diagonal  aatrlz  ■hose 
eigenvalues  are  the  aagnltudet  of  the  eigen* 
values  of  T.  This  aatrlz  represents  a  trace 
T  on  t. 


das  operator 


2.2.  Recursive  Definition  of  the 
eoefflclents  v  ,t  T 

For  1  •  I.2.3.**',  *0  defint 


ARHA 


tj.  and  the  I, 's  in  (3.2) 

ore  the  coefficients  generated  by  (2.3)  and 
(2.4)  for  T.  Of  course. 


|a,|s  I,,  1  -  1,2, 


(3.3) 


e,  ■  tr  S^/1 


(2.3) 


For  each  n 


(1) 


let  q 


(1) 


be  the  sua  of 


■here 


$,  ■  T  and  S,*,  •  TlSi'a,!). 


*1*1 


r*i‘ 


(2.4) 


Alto. 


the  absolute  values  of  all  eigenvalues  of  T 
eicept  those  In  l^)  .  Hence 

qj^’sl,  .  1,k.l.2.....  (3.4) 

The  hay  Identity  here  It 


b,  •  a  and  b^  •  C(S,.,-a,.,I)8.  1  a  2.  (2.S) 

These  a^'s  and  b^'s  are  the  coefficients 

that  ue  use  In  (1.3).  If  the  state  space  C 
hat  finite  dlaentlon  n,  then  a^  ■  b^  ■  0  for 

1  >  n.  In  this  cate,  the  a^'s  are  the 

negatives  of  the  coefficients  In  the  charac¬ 
teristic  polynoalal  of  T,  and  (2.3)-(2.S) 
constitute  a  uell  knoun  algorltha  forcoaputlng 
the  ARHA  coefficients.  The  algorltha  here 
then  uould  appear  to  be  the  natural  uay  to 
define  the  ARHA  coefficients  for  the  Infinite 
dlaenslenal  cate.  Houevar,  that  the  output 


I  •  (Id)  I  •  1.2.—.  (3.5) 

k«,  •  *  k*l  “ 

Froa  (3.2).  (3.4)  and  (3.5).  «c  have 

■id  ‘TO  ‘l*  ’  ■  I.3.*-.  (3.6) 

■hich  yields 

*1  •  ®  O  ’) 

for  any  r  >  0.  Then  the  theorea  folleas 
froa  (3.3)  and  (3.2). 


Our  nuatrlcal  tiptrltnct  nUH  utapUs 
tugttits  tkat  tkt  coafflcttnta  k|  alto  Oocay 

fattar  that  any  aaponanttal.  undar  our 
kypothatat  on  T.  t  and  C.  So  far,  kouauar,  ho 
can  prove  tkit  only  wktn  tka  aloanvoctort  of  T 
for  an  ortkenoraal  katit  for  E.  In  tkit  cato. 
an  arguoent  tiollar  to  tka  proof  of  Tkooras  1 
tkowt  tkat  II  S^llj  dacayt  fattar  tkan  any 

anpenantlal.  Ma  kopa  to  oatand  tkit  ratuU 
for  any  T  of  traca  data. 


*1n  "  ^In  “  ®'  I  •  •♦It  •♦i.**'.  (4.9) 

Oacall  tkat,  bafora  now.  It  kat  not  bean 
thown  tkat  tka  y(t)  of  (1.2)  tatltflas  (1.3). 
If.  for  tooa  t^.  wa  choota  z„(tg)  tuck  tkat 

11o  *„(t^)  ■  aCtj).  (4.10) 

tkan  (4.1)>(4.3)  along  witk  (4.C)>(4.7)  aniura 
tkat 


Mkan  tka  tpaetral  radiut  of  T  It  latt  tkan 
1,  wo  can  tkow  at  laatt  tkat  b|  approackat  taro. 

Tkit  appllat  to  tka  coapantator  of  tka  oxaapla 
In  Section  S. 


IV.  toprotlaatlon 

Ma  attuaa  tkat  wa  kava  taguancet  of  finite 
rank  oparatort  T||,  0^,  tuck  that 

l|T„-T||,-0.  (4.1) 

II  »„-l|l  ♦  0.  (4.2) 

II  veil  *  0.  (4.3) 

at  n  ■•  ».  for  aack  n,  wa  cowpute  ARHA 
coafflclentt  a^^  and  b|||  according  to  (2.3)-(2.S). 

Than  (4.1)-(4>3)  guarantee  tkat,  for  aack  1. 

1^  t,„  -  a,  (4.4) 

and 

1Ubin-bi.  (4.S) 

In  appllcatlont.  (4.2)  and  (4.3)  tlwott 
tlwayt  hold  bactuta  •  and  C  have  finite  rank. 
Mowaver.  (4.1)  It  ttrongar  than  tka  conworgonca 
obtained  with  atny  approxlwatlon  ackawat,  and 
wa  kopa  to  be  able  to  weaken  tkit  condition. 

Our  approxlwatlon  of  tka  opan>1oep  tyttaw  In 
tka  exawple  of  the  next  taction  llluttratat 
tka  Iwportant  clttt  of  appllcatlont  wkara  (4.1) 
koldf  bactuta  the  e1|anvectort  of  T  forw  an 
ortkonorwal  batit  for  E  and  tka  approxlwatlon 
contitti  of  projecting  onto  the  algantpacat  of 
t  (wodal  approxlwatlon). 

Alto.  If  the  trantitlon  operator  It 
T  ■  Tg  ♦  F.  where  the  alganvectora  of  forw 

an  ortkonorwal  bttit  and  F  kat  finite  rank, 
tkan  approxlwatlon  by  projecting  onto  the  olgan* 
tpacet  of  Tg  will  yield  (4.1).  Tkit  It  the 

cate  In  our  approxlwatlon  of  tka  cowpentator 
In  tka  oxawpla. 


■  Vn<*>  ♦  V'*** 

nenslonal 

(4.d) 

»n<*>  •  Ve<t). 

(4.7) 

• 

(4.0) 

wkara 


11w  y  (t)  -  y(t).  tat.  (4.11) 

n»» 


If 

y(t)  •  y„(t)  •  u{t)  •  0,  t  <  t^,  (4.12) 

than  (4.1),  (4.2)  and  (4.11)  Iwply  that  aack 
tarn  In  (4.8)  convargat  to  each  tarn  In  (1.3). 
Nltkout  (4.10)  and  (4.12).  wa  justify  (1.3) 

In  a  tlwllar  way,  but  only  at  t  •>  >.  Tkt 
decay  rates  for  a^  and  b^  are  tkan  useful. 

V.  Aopllcatlen  to  Flexible  Structures 

S. 1 .  Abstract  Structure  Model 

An  Iwportant  class  of  applications  for 
Infinite  dlwtnslontl  AANA  wodals  Is  the  control  of 
flexible  structures  whose  linear  distributed 
wodal  hat  the  forw  of  tka  differential  aqua¬ 
tion 


V  ♦  "o*  *  *0*  •  V* 

where  the  stiffness  operator  A„  usually 

contains  partial  diffarantlal  operators  for 
the  elastic  conponentt  of  the  structure.  In 
applications  A^  Is  salfadjoint  with  cowpact 

resolvent  and  salfadjoint  and  coercive. 

For  details  of  tka  first  order  forw  of  the 
aquation,  sea  [3].  Ha  only  note  tkat.  with 
roasonabla  conditions  on  tka  dawplng  operator 
0^,  wa  can  taka  tka  state  vector 

s(t)  -  (x(t).  xlt)) 

and  write  tka  solution  to  (S.1)  at  (1.1)  for 
pfecawita  constant  control  u(t).  Of  course, 
wa  scale  the  tlwe  variable  so  tkat  the  Input/ 
sawpllng  Interval  Is  1.  For  tka  observability 
condition  tkat  we  sssuwad  In  Section  1,  It  It 
sufficient  tkat  no  natural  frequency  of  the 
wndswpad  structure  be  a  wultipla  of  the  sawpllng 
frequency. 

A  cowwon  dawplng  wodal  for  flexible 
structures  It  wodal  dawplng  that  provides  the 
sane  dawplng  ratio  for  aack  woda.  This  weans 
(sea  [1])  tkat 

Oo  ‘  *0*5 

In  (S.1),  with  c^  a  positive  real  nuwbar.  It 

can  be  shown  that  •  at  lease  for  structures 
whose  flexible  cowponents  are  beans,  plates, 
strings  or  wewbranes  •  the  dawplng  operator  In 
(S.2)  causes  the  operator  T  In  (1.l)  to  be 
trace  clast. 


Me  should  note,  however,  tkat  linear 


0-  •  e  A 
0  0  0 


(ttt  [3]) 

(5.3) 


products  a  non-trtct  clast  T  fetcauta  a  tubta* 
putnct  of  tha  alganvalutt  of  T  convargat  to 

For  (5.1)  with  althar  (5.2)  or  (5.3), 
tha  alganvactors,  or  *0500 ,  of  tha  opan>1oop 
tpttaa  fora  an  orthonoraal  basis  for  tha  stata 
spaca,  so  that  aodal  approilaatlon  ylalds  (4.1)> 

(4.3) .  Mt  usad  aodal  approilaatlon  In  tha 
folloaing  aiaaple,  but  we  ora  particularly 
Intarattad  In  being  able  to  aeaktn  (4.1)  to 
justify  aora  general  finite  aloaant  approilaa* 
tions  for  approxiaating  ARHA  aodals  of  flexible 
structures. 

5.2.  Exaaole 

Tha  structure  In  this  axaapla  consists  of 
an  Eu1ar>Bernou1 1 1  beaa  cantlllvarad  to  a  rigid 
disc  which  Is  free  to  rotate  about  Its  fixed 
canter.  In*p1ana  aotlon  It  aodalad.  Including 
linear  transverse  vibratipns  of  tha  beaa.  The 
rotation  of  tha  disc  gives  a  rigid-body  aoda. 

An  actuator  applies  a  control  torgue  to  the 
disc,  and  a  tensor  aaasurat  tha  rotation  of 
tha  disc. 

For  this  problaa,  we  utod  too  levels  of 
approilaatlon.  First,  we  used  a  finite  elaaant 
approilaatlon  to  coaputa  tha  natural  aoda 
thapas  and  fraquanclas  of  tha  coaposita  tystea. 
Than  we  used  tha  tyttaa  aodas  at  tha  basis 
vectors  In  tha  approxlaatlont  of  tha  ARHA 
aodal  and  tha  optlaal  control  problaa.  For 
coaputing  both  ARHA  coafficlantt  and  control 
and  attlaator  gains,  wa  found  that  tha  first 
25  tyttaa  aodas  wart  aatlly  sufficient  for 
eonvargance.  Me  aada  sure  that  our  original 
finite  elaaant  approilaatlon  had  convergad  for 
tha  first  25  aodas. 

Tha  lowest  natural  freguancy  of  tha 
structure  was  4.f  rad/tac.  Me  used  the  aodal 
daaping  aodal  In  (5.2)  with  both  c.  ■  .02 
(15  critical  daaping)  and  c^  ■  .1  ”  (SI  critical 
caaping). 

No  coaputad  the  open-loop  ARHA  coefficients 
at  ditcuttad  In  Section  4,  with  n  being  the 
nuabar  of  tystea  aodas  Included  In  tha  approxl- 
aatlon.  Table  1  shows  tha  converged  (recall 

(4.4)  and  (4.5))  values  of  tha  ARHA  coeffi¬ 
cients,  roundad  to  two  significant  digits. 


1 

.43 

Ea1 

.29 

E-2 

2 

-.07 

E«1 

-.11 

E-1 

3 

.11 

E+1 

.21 

E-1 

4 

a 

-.02 

E«1 

.25 

E-1 

a 

10 

.30 

E*1 

-.16 

E-1 

11 

a 

.39 

Eel 

.17 

E-1 

a 

20 

.32 

1*0 

.19 

E-2 

21 

a 

.04 

E40 

-.35 

E-2 

30 

-.40 

E-1 

.14 

E-3 

31 

a 

.32 

E-1 

-.90 

E-4 

a 

40 

-.39 

E-S 

-.11 

E-7 

41 

a 

.70 

E-6 

-.22 

E-0 

SO 

-.11 

E-12 

.36 

E-1S 

51 

.14 

E-13 

-.36 

E-16 

Table  1 


1 

.26 

t4l 

.29 

E-2 

2 

-.22 

Efl 

-.62 

E-2 

3 

-.51 

E^1 

.30 

E-2 

4 

• 

.16 

Eel 

.34 

E-3 

lO' 

-.00 

E-1 

-.21 

E-2 

11 

a 

.24 

EaO 

-.37 

E-3 

to 

-.15 

E-3 

.44 

E-6 

21 

a 

.21 

E-4 

-.50 

E-1 

a 

30 

-.51 

E-14 

.16 

E-16 

31 

.47 

E-16 

-.07 

E-10 

Table  2 


il  106  Control 


For  (1.1),  we  use  tha  perforaonce  Index 

J-I  tll*(t)ll*  (5.4) 

t*0 

where  a  >  0.  The  optlaal  control  law  than  It 
u(t)  •  -K  i(t),  (5. 5) 


«  •  [’♦•JW.]*’  •JfT,  (5. 5) 

and  the  nonnogatlvo  tolfadjoint  oparstor  F 
satitflat  the  Infinite  diaanslenal  Riecatl 
agnation 

f  •  I  ♦  T*(R.»,tl4t;p»J-\p)  T^,  (1.7) 


■  i.'T^  • 
M  ,V 


■her*  T  ■  T/a  and  ■  B/b.  A*  usual,  the 

S  tt 

••shift  In  this  probItB  ansurts  that  tha 
spactral  radius  af  tha  c1osad*1aep  aparatar 
T>BK  Mill  ba  no  graatar  than  a.  In  this 
axaapla.  va  usad  b  ■  .B. 

Tha  stata  astlMter  far  (1.1)  has  tka 


i{t*1)  ■  T  i(t)  ♦  Bu(t)  ♦  «(y(t)-C  i(t)).  (S.B) 
To  gat  tha  gain  6,  Ma  solva  tha  optlaal 


ylalds 

•hara  P  satisflas 


'a  ■  •’ 


I 

-.14  Ea1 

.31 

E«3 

tka 

2 

-.44  EfO 

-.S2 

Ea3 

3 

.17  Ea1 

.17 

E«3 

4 

•.9S  E40 

.19 

Ea9 

($.0) 

10 

.78  E-1 

-.17 

$♦1 

the 

11 

-.81  E-1 

.24 

Ev2 

s 

20 

.1$  E-3 

-.41 

E-1 

($.9) 

21 

• 

-.19  E-4 

.94 

E-2 

• 

30 

.$$  E-14 

-.16 

E-11 

($.10) 

31 

-.22  E-1$ 

-.20 

E-14 

and  C_  ■  C/b. 

B 


In  stata  spaca  fora,  tha  eaapansator  Is 
t(t«1)  -  (T-fiC)2(t)  *  Bu(t)  a  t  y(t)  (S.11) 


Tabla  3 


ifaraneas 


Mara,  tha  a*th1ft  arsuras  that  tha  spactral 
radius  of  c1otad>1oop  astlaator  operator 
T*6C  Is  no  graatar  than  a.  Again,  Me  usad 

B  «  .9. 

Mom,  for  (S.11)<(S-12).  Ma  think  of  y(t) 
at  the  Input  to  the  eaapansator  and  u(t)  at 
the  eutput.  Applying  tha  aethad  of  Section  z, 
Ma  obtain  tha  AANA  raprasantatian 

a(t)  ■  I  [a5u(t-1)  ♦  b?  y(t-1)].  (S.19) 

1»1  ' 

To  coaputa  tha  Oi't  and  bf's  far  the 
asaapla,  mo  solved  a  iaguance  'of  finite 
diaantlonal  IQ6  probleat  corresponding  to  tha 
toquanca  of  aodal  appreilaatlons  to  aur 
flasibla  structure.  (See  [A]  far  details.) 
This  yields  a  sequence  of  finite  diaantlonal 
ceapansators  of  the  fora  ($.11)>(S*12).  Far 
each  such  eaapansator,  mo  caaputad  the 
corresponding  ARNA  aodal  of  tha  fora  ($.13) 

Mlth  caofficlantt  a^^  and  b|^,  Mhara  n  It  the 

nuabar  of  aodat  In  tha  systea  approslaatlen. 
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1  •  1.2,- 


($.14) 


(S.IS) 


coefficient  It  c.  ■  .1.  Table  3  gives  tha 
ARNA  coefficlantt  for  ($.13). 
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